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— Introduction . The information about energy spectrum of sea 
level muons at high energies beyond magnetic spectrographs can 
be obtained from the underground intensity measurements if 
the fluctuation problems are solved. In the present paper we 
recalculate the correction factor R for the range fluctuations 
of high energy muons by analytical method of Zatsepin et. al 1 ) 
where most probable energy loss parameter are used. It is 
shown that by using the R at great depth together with the 
slope, A, of the vertical depth-intensity (D-I) curve in the 
form of exp(-t/A), the spectral index, y, in the power law 
energy spectrum of muons at sea level can be easily obtained. 

2. Formulation of the correction factor R. The R is defined 
as R = lQ(t)/I(t) , Iq ( t) is the intensity without fluctu- 
ations and I (t) is the intensity , taking fluctuations into ac- 
count. In order to find I (t) the following diffusion equation 
has to be solved: 

3I/3t - B(E )3I/3E = f'w(E ,v)[I(E /l-v,t)-I(E ,t)]dv (1) 

” " J o ^ U V 

W(E ,v) = (N/A) $ + N0 . t 

j H £5 iN 

(N: Avogadros number. A: atomic mass number), where $ and 

B 

$N are the differential cross section of bremsstrahlung and 
nuclear interactions, respectively. To get the solution of 
eq. (1) , we make two assumptions and introduce three dimension- 
less variables as follows: 

(i) B(E y ) = a + b E , (ii) f(v)/v = W/(b B +b N ) 

X = V' ' = V/ a - b * (b B +b N )/b p- 
, where a is ionization loss and b^, and b^ are energy loss 
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parameters due to pair creation, bremsstrahlung and nuclear 
interaction, respectively. Then eq. (1) becomes 

dj/dx - (e + 1 ) 3 1/3 e = b f f (v) /v [I ( e/l-v , x) -I ( e , x) ] dv (2) 

J 0 

In this notation the average energy loss may be written as 
-de/dx = 1 + (1 + b)e (3) 

Iq(e/x) is given by combining eq. (3) with the boundary con- 
dition l(e,0) = Be~ Y as 

I Q (e,x)= B-exp [-y (l+b)x] [e+ (l-exp(- (l+b)x))/ (1+b) ] _Y (4) 

By analogy of eq. (4) we take the solution of eq. (2) in the form 

I (e,x)=B* exp[-y (1+xb) x] [e + (1-exp (- (1+Kb) x))/ (1+Kb) ] _Y 

x exp [<(> ( e , x) ] (5) 

From eqs. (4) and (5) we get R for e = 0 in the following way 

R(0,x) = [exp (b(x-K)x)-(l+b)/ (1+Kb)* [exp ( (1+Kb) x) -1)/ 

[exp ( (1+b) x) -1] ] +Y exp [-<J> (0,x) ] (6) 

2 

At a great depth(t> 4000 hg/cm ) and under the assumption 
4> (0, t+°°)=0, we get such a simple form as 

R [ (b p + b f)/( b p + Kt> f ) ] Y exp[-Yb f (1-x) t] (7) 

, where b, = b^ + b„. 

f B N 

3. Derivation of spectral index y. It is assumed that the 
D-I curves are expressed by a unique exponential law of the 
type I (t) = C*exp(-t/A) (8). Then combining eqs. (4), (7) 

and (8) we can obtain a relationship between y and A : 

Y = [A (bp + x b f )] _1 (9) 

Here the constant x is determined in such a way that i)>(e,x) 
approaches zero as e+°°. By substituting eq. (5) into eq.(2) 

X is given as 

/I 

X = (1/ Y ) ) f (v)/v[l-(l-v) Y ]dv (10) 
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For 4> , $ we take the following formulae 

D N 

$ B =0.95a(2r 0 m e /m vi ) 2 Z(Z+4) (4/3+v 2 - (4/3) v) (l/v)ln(P/Q) (11) 

P= (2/3) k (m /m ) Z -2/3 , 

M “ 

Q=(k/e-/2) (m 2 /m E) (v/(l-v))Z _1/3 + 1 2) 
m “ y 

{ ( v-l)-|- [l-v + y(n-^)] 

r 1+ i 3) (i 2 ) 

x i n — — L — zShiJ l 

Ev /Ev . AN 
/f V A 2M/ 

(a^ N = 125 pb and A 2 = 0.4 Gev 2 ) 

We calculated the values of x as a function of y for 
standard (Z=11,A=22) and K.G.F (Z=12 . 93 ,A=26 . 12 ) rocks at E =1, 
10,100 TeV. The results are shown in Fig.l. We find from this 
figure that the x and y has nearly a linear relation such that 

X = m 1 Y + m 2 (13) 

In the case of E =10 TeV we have m, =-0.1085, m„=0.9943 for . 

y 12 

S rock and m^-0.1092, m 2 =0.9935 for K.G.F rock. Substituting 
eq. (13) into eq. (9) we get the equation for y . 


m, bj-y + (b + m 0 b f )y -100/A = 0 (1< 

1 r P z r _6 -2 

(b ,b f are measured in units of 10 cm g and A in units of 

P 1 lft 4 -2, 

10 gem ) 


For b , we use the expression taken from Bugaevs booK 
P 

b =6.01*10 -8 Z(Z+1)/A[0.971n(p 1 /Q 1 ) + 2 . 15 ] (g -1 cm 2 ) 


p-^isiz 


-1 / 3 , Q 1 =183Z _1 ^ 3 m 2 /(2E m ) + 1 

' 1 y ye 


In the following we apply this equation together with the A 

of D-I curves for S rock and K.G.F rock to get y If we take 

4-2 -2 • 

A=9. 868*10 gem (8000 ^t ^9000 hgem ) for S rock, which is 

derived from eq. (2) in ref. (5) , and A=9. 00*10^ gem 2 (t=10^ 
hgem 2 ) for K.G.F rock^ 
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The values of y thus obtained is as follows: 

Y = 2.70 (E 3 0 TeV) for S rock and 
y = 2.62 (E*»50 TeV) for K.G.F rock. 

The former is consistent with the value of 2.71 by Bergamasco 
et al. 5 ^ and the latter is also consistent with that of 2.6 
by Miyake et al."^, where fluctuations are treated by Monte 
Carlo Method. 

4. Conclusions. By using the approximate R at great depth, y 
can be easily obtained if the D-I curve has the form I (t) = 

C • exp (-t/A) . It is found that there is a discrepancy of y 
between S rock and K.G.F rock in the same energy region. 
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